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ABSTRACT 


Specification of parameters defining the mean, as well as the 
turbulent properties, of the over water surface layer is examined. A 
baroclinic boundary layer model is considered using Cardone's modifi- 
cations for the marine environment. 

Data fields of pressure and temperature were selected as input, 
and original ship observations were utilized to verify the predicted 
surface parameters. The effects of stability changes, isobaric curva- 
ture, thermal wind, and the value of the von Kármán constant in the 
model were examined with results shown. 

Iteration criteria used in the model were evaluated with respect 
to operational application. Winds calculated by the model at the height 
of 19.5 meters were incorporated into the current FNWC wave height 
program, and the results compared with observed wave data. Turbulent 
parameters describing the properties of the index of refraction CA were 


also computed and evaluated with respect to laser propagation criteria. 
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I. INTRODUCTION 


Parameterization of the boundary layer over the sea surface has 
been receiving increased attention in recent years. Fortunately much 
of the theory developed for low level turbulence over a land surface is 
also applicable over the-oceans. But the air-water interface does pre- 
sent some unique problems. For instance a fluid surface, rather than a 
fixed one, results in mobile surface roughness elements. Because of 
the presence of a fluid surface, any parameterization requires a knowl- 
edge of the processes which transfer heat and momentum energy across 
the boundary and ultimately influence such variables as wave generation 
and index of refraction. 

The purpose of this study was to examine and evaluate several 
features of a recently formulated marine boundary layer model which was 
designed for operational use. The model is a two-layer baroclinic model 
for the marine boundary layer and the formulation examined was that 
developed by Cardone (1969). Cardone based his formulation on an 
earlier two-layer model due to Blackadar (1965a). Features of interest 
in the model are related to the physics of the planetary boundary layer 
and include the effect of atmospheric stability, baroclinicity, coriolis 
acceleration and changes in the surface roughness parameter on the 


pressure driven surface wind. 
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The scope of the study was to examine the model in the context of 
its operational usage. As such, it was a case study in which the data 
represented pressure and temperature observations from sequential 
12-hour periods when a large cyclonic pressure system dominated the 
Central North Pacific Ocean. This particular synoptic situation provided 
the necessary relations between the scalar fields to examine features 
associated with baroclinicity, stability and wind speed categories. 

Because these data represent scalar fields which had to be de- 
fined by some analysis scheme before being used as input to the two- 
layer model, the differences which arise due to procedures in the external 
field specifications are considered. The effect of isobaric curvature on 
the surface wind fields is examined and this effect is compared to the 
changes associated with the physics in the model. The differences 
which may arise solely because the fields are defined by an objective 
versus a subjective analvsis are considered. 

Finally, expressions describing the boundary layer are, by their 
very nature, empirical. As such, they contain empirical constants and 
the von Kármán constant (k) is one whose value has recently been ques- 
tioned. The importance of a suggested change from k=0.40 to 0.35 was 


considered. 
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II, ASPECTS OF A PLANETARY BOUNDARY LAYER 


A. SURFACE BOUNDARY LAYER 
1% Basic Structure of the Boundary Layer 

Recent advances in the study of atmospheric turbulence 
have led to a greater confidence in describing the distribution of wind 
and temperature in the atmospheric boundary layer over a solid suríace. 
This progress has arisen largely through the application of the Monin- 
Obukhov similarity theory which has been supported by recent observa- 
tional studies. Cardone's contribution was primarily to extend 
Blackadar's two-layer neutral baroclinic boundary layer model over fixed 
terrain, and then to apply this model to the boundary layer over a fluid 
surface. The major difference between these boundaries is that the sea 
surface changes its character as the wind blows over it. 

In general, Blackadar's model treats the boundary layer as 
being formed by both turbulent friction and the coriolis acceleration. All 
processes are considered to be three dimensional, and the layer is char- 
acterized by a density stratification. It is conveniently separated into 
the three regions depicted in Figure 1. In this delineation the important 
consideration is the specification of Kj, which is the turbulent transfer 


coefficient for momentum. 
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Free Atmosphere 


€ 600 m 
۱ 
۸ | Ekman Layer 
2 ۱ 
h «50m 
A Surface Layer 
L^ Viscous Sublayer 
0 TT a əəə <i cm 
ce 


Fig. 1. Diagram of the planetary boundary layer 


a. Viscous Sublayer (less than 1 cm) 
In this near surface layer the stress is supported by 
viscosity (aerodynamically smooth flow), and the mean velocity 
increases linearly according to 


— a * (1) 
U: 3 ۷ 





where v is the kinematic viscosity and U, is the friction velocity. 
b. Surface Laver (up to 50 meters) 
This is the constant stress layer in which the coriolis 
parameter can be neglected and Kn increases linearly with height. For 


neutral stratification the logarithmic law is valid 


16 





U 


هذا 
[z k log / 2 )2(‏ 


where k=0.4 is the von Kármán constant. 

e Ekman Layer (50-600 meters) 

This layer is characterized by the stress decreasing 
with height and К constant throughout the layer. 

In formulating the boundary layer model, Cardone (1969) 
considered the validity of the "log law" for the surface layer over water, 
which has been experimentally documented down to near surface limits 
(Phillips, 1966; Paulson, 1967). He also assumed that the process of 
wave generation did not seriously affect the surface layer since under 
neutral and steady conditions a logarithmic shear zone should be estab- 
lished. Accepting the log profile as being valid, Cardone was able to 
apply much of the surface layer turbulence theory developed over land 
to over-water conditions. However Cardone did take into account the 
action of the waves when specifying the effective roughness parameter 
for the sea surface. 

2. Similarity Theory for the Surface Layer 

The Monin-Obukhov "similarity theory" has led to consider- 
able progress in the specification of the distribution of wind and tem- 
perature in the boundary layer. The basis of this theory is that near 
the ground there exist velocity (U,), length (L), and temperature (9,) 
scales which are essentially invariant with height. When the principle 


variables such as temperature (T), wind (U), and height (z) are expressed 
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as nondimensional fractions of these quantities, a series of nondimen- 
sional equations can be formed which are of general validity in the sur- 
face boundary layer (Lumley and Panofsky, 1964). Recent analyses of 
the oceanic wind data (e.g. Phelps and Pond, 1971) have indicated that 
the Monin-Obukhov theory is valid for over-water application. These 
results support Cardone's inclusion of the Monin-Obukhov stability 
criteria in the boundary layer model. 

Similarity theory predicts that a universal relation should 


exist for the nondimensional wind shear 


2 9 
Ux SE (3) 


and the nondimensional temperature gradient 


2 bo 
Ox IZ | 0 (4) 


where J, and gy are functions of the dimensionless ratio z/L. Here 


9 is the potential temperature, 9. is the scaling temperature defined by 


1 H 


Cp Pa P‏ يپ 


Ox 
and L is the Lettau-Monin-Obukhov stability length defined by 


3 — 
_ -U& Cp PaT 


The heat flux (H) is not easily measured, so a modified stability length 


was defined by Panofsky (1963) 


Ky UT ou ez 


E 00 9۹9۹ ۷ ° 


(7) 
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The assumption that L' is independent of height implies that ар is con- 


stant, which in turn implies similarity between wind and temperature 








profiles. 
If L'is a valid scaling length, similarity theory predicts 
that 
kz UM | 
U : öz = an (1221) (8) 
Ox öz ah : )9( 


Integration of these formulae (Panofsky, 1963) yields 





Ue — ۲ 254 0ت2‎ (10) 
9 
9289 + 5 E 22 د‎ | , )11( 
where the relation between W and 2 is 
z/L 1-g (E) 
V (z/L') - mee SE Wet ur. (12) 
O 


Paulson (1970), using the KEYPS formula, integrated (12) by parts to 


obtain the following expression for W 


٦ 


۱٣۱۰۰٥۳2000 >‏ 1214 - لہا 





E 2 
+ 2 tan “A - — + In Ei (13) 


Within the scope of an operational application, these 
expressions for the surface layer would be used to determine the wind 
speed given Ux , or to determine Ux given the wind speed at a particular 
level. For this reason the partial derivatives in (7) are replaced by 


expressions involving height differences and air-sea temperature 
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differences. This is achieved by substituting the following expression 
for M , Obtained by combining (9) and (11) 


Z 
do. B, (T) (8, - 8.) 
öz ” z Lin(73/z9)- Q Za 7L] ' . 


where z. is the height at which 94 is measured (chosen as 10 meters by 


a 
Cardone). 


Substituting (14) and (8) into (7) yields 


| vé q [1n 22) - ф(2а/т,)| 
L' = ESG тә e + (15) 


Equation (10) is then solved for Ux 


k Um 
۸ وي‎ )16( 
A m) 


where 27 is the height at which the wind speed is measured (chosen as 
19.5 meters by Cardone). Equations (15) and (16) can then be solved 
for Ux given the wind speed at some level (U,) and the air-sea tempera- 


ture difference. 


3; Consideration of the Drag Coefficient ( Cz) and Roughness 
Parameter (70) 


The surface roughness parameter and the influence of thermal 
stratification have received considerable attention in the definition of 
the relation between the nondimensional drag coefficient (C4) and the 
wind speed at a certain height (usually 10 meters). The need for such a 


relation arises from the "bulk aerodynamic" formula for estimating 
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momentum flux. This equation, deduced from dimensional analysis, has 


the form 


r= Pp, Uk = p. CG, U2 , (17) 


where T is the surface stress. From this expression the drag coefficient 


is defined as 





o E (18) 


Rearranging and squaring (10) yields the following expression for the 


drag coefficient at a level z in terms of 2 and stability 


5 2 
"7 [m 7/zo - w(C/L) 2 
which was the form considered by Cardone. 
Observations have indicated that over the sea, z, is depend- 


O 
enton U,. This dependence, however, changes as the wind speed 
changes. For the case of light winds (aerodynamically smooth flow), 
the drag coefficient 2 decreases with increasing wind speed (Kraus, 
1966) and zo has the form 

۷ 


Soo ии (20) 


For aerodynamically rough flow, 5 increases vvith increasing vind 
speed according to 


toT کي‎ 22) 


where (a) is a constant = 0.035 (Charnock, 1955). A summary of field 
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observations indicates that the change between rough and smooth flow 
occurs at wind speeds of 5-6 meters/second. 

Cardone proposed the following simple expression for z 
covering a wide range of wind speeds (both rough and smooth flow) 


zı 


O U, 





+ UE “5. (22) 


The constants vvere chosen so that ce was a minimum fora value of 
6 meters/second at the 10 meter height under neutral conditions, and 
above this speed the relation approximated Charnock's results. For Ux 


in meters/second and z. in meters, the resulting expression was 


0 


3 


9 = - 4 
Zo = 6.84 x 10 ur + 4,28 x 10 us - 4.43x 10, (23) 


Other representations of C, versus U, have been published. 
A few of these appear in Figure 2 along with Cardone's approximation. 
It is clear from the scatter of curves in Figure 2 that the specification 
“of the drag coefficient was a primary decision in Cardone's formulation. 
This choice of us Versus سلا‎ will be examined during the discussion on 
results from this investigation. 

With the equations for stability length (15), friction velocity 
(16), and surface roughness parameter (23), the surface layer wind dis- 
tribution is completely described within the framework of the similarity 
theory. These three equations can be solved simultaneously for Ux 
using iteration techniques with the external parameters described below 


as input. 
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Fig. 2. Various suggested forms of the drag 
coefficient (C,) versus windspeed (ue) 


at 10 meters 
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B. PLANETARY BOUNDARY LAYER 

7 Background 

In developing a model for the entire boundary layer, con- 

sideration has to be given to the possible external parameters. In 
Cardone's formulation these гара were the geostrophic wind (G), 
the air-sea temperature difference (0. 5 the coriolis parameter (f), 
and the thermal vind (r). He considered previously suggested empirical 
relations between the geostrophic drag coefficient (U,/G) and the sur- 
face Rossby number (Lettau, 1959) 


٦ G 
NE EZ 





; (24) 


where G is the surface geostrophic wind speed. Empirical relations 
have also been suggested between the cross isobar flow angle Y and 
Re 


O 
Blackadar (1965a) developed a two-layer boundary layer 
model, consisting of a Prandtl surface layer and an Ekman layer, whieh 
was verified on the basis of the above relations. 
Cardone's formulation was closely related to Blackadar's 
work. In particular, the height of the surface layer (h) was specified 


explicitly in terms of the external parameters 


Bo G 


h = f 


(25) 


: -4 
where By was a nondimensional constant equal to 3.0x 10 as deter- 


mined by Blackadar. 
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By applying proper boundary conditions and requiring that 
continuity of wind and wind shear exist at level h, the following solu- 
tions for U&/G and YW, (along with the solution for 2.) completely 
define the surface boundary layer wind distribution under neutral 


conditions. 





U 3/2 
* 
E 55 ə ve | (26) 


Us kY2 sin (n/4- Wo) 


“ur E O (27) 
G 


In B se 





Implicit in such an approach is that the value of К at the internal 


boundary between the surface layer and the Ekman layer is given by 


KU: G 


5 (28) 
11 f 


2, Cardone's Extension of Blackadar's Two-Layer Model 


Ihe success of Blackadar's two-layer model in predicting 
the essential characteristics of the neutral surface boundary layer, along 
with its relative ease of application, suggested its extension to the non- 
neutral marine boundary layer. Cardone extended the diabatic model to 


over-water flow by regarding z, as a function of Ux. 


O 

The eddy viscosity in the constant stress layer for the non- 
neutral case was expressed by Cardone in terms of the nondimensional 
wind shear as 


8 kU«z 
En = . (29) 


ДИ) 
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Since K,, is constant with height above h, its value in the Ekman layer 


is then 
kU, B.G 
۱ rN; 
The wind speed at h was given by 
U 
Up = m 700 9 ٤ (31) 


Internal boundary conditions applied at level h appear in Figure 3 and 
follow from the fact that 5 always makes an angle of ru. to 
the ageostrophic wind (W) and that ou is parallel to U in the 


constant stress layer. 





G 


سے 


Fig. 3. Wind vector diagram applicable at level h. 


By requiring continuity of wind shear at h, as wellas directional 
continuity, Cardone was able to reduce the problem (analogous to (26) 


Sud (27)) to 








Ux | 2 1/2 

E Ë k Bosin WB, 2221 (32) 
Ја  КлЈ2 sin (1/4-U,) ə) 
G in Bo Ro - Y (h/L;) 
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U 
These equations are solved simultaneously for Er and Vo from the 


external parameters G, f, and (8.-9 ). In addition, a dimensionless 


S 


stability parameter was defined as 


[NM 
k ET (34) 





which represented the ratio of the surface layer h to the stability 
length L' (Blackadar and Ching, 1965). 

The relationships determined by Cardone for U«/G and УД 
as functions of Ro and L, appear in Figures 4 and 5. For all surface 
Rossby numbers, U,/G increases and Wo decreases with decreasing 
stability. In general, stable conditions have a more pronounced effect 
on the suríace layer wind characteristics than unstable conditions. In 
comparing the results from the model with observations, Cardone noted 
that Ux/G increased by roughly 30% in unstable conditions, with no 
marked tendency to decrease at great instability, and decreased roughly 
70% in stable conditions. These features in the observed results will 
be discussed later. 

The effects of baroclinicity were also included by assuming 
that the geostrophic wind was a linear function of height (Blackadar, 
1965b). Blackadar (1965a) introduced a parameter, p, which was defined 
at level h by 


UL B¿G 


hyr:‏ - ا 
P = 190/32), ^ fA, (I) b cree | "E‏ 
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Fig. m Variation of the geostrophic drag coefficient 
with surface Rossby number and dimensionless 
stability length (L,) (Cardone, 1969) 





logo Ro 


Fig. 5. Variation of surf&ce cross isobar angle vith 
surface Rossby number and stability (Cardone, 
1969) 
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The requirement of continuity of wind shear for the baroclinic case took 


the form 


U, = | 99), Al (36) 


Figure 6 illustrates the geometrical implication of (36). From the Ekman 


layer solution of the ageostrophic wind, the following equation can be 


E 
P öz 


where p' is a dimensionless parameter defined by 


BG E (B/1) 3 
p! 9 7 25086 E 5077 v e. (38) 


obtained 


= N2 p | w, | : (37) 
h 








Fig. 6. Geometrical relationships resulting from the 
requirement of the continuity of wind shear at 
the internal boundary level h (Blackadar, 1965b). 


Next the thermal wind vector was replaced by its dimensionless 


magnitude 


- i |] 
ns.” öz 


0 








(39) 
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ð G 
where the direction of Sr relative to the surface geostrophic wind 
was implied by angle n in Figure 6. Still another dimensionless 


quantity was defined 
i = h ' 
r E b LL”) (40) 
and (35) was then rewritten as 


سے 
Pi öz‏ 


By applying the law of cosines to each of the two smaller triangles in 


eu er (41) 





Figure 6 that have M for one side, the following relations were determined 


2 

RES | лг, | 5—: (42) 

s” z (132p"32p”) (43) 
2 2 

a = (1+r'“-2r' cosn ). (44) 


n is the angle between the geostrophic wind and the thermal wind (see 


Figure 6). Then the law of sines was applied to yield 


sin d - 9 ع/‎ (45) 
sin Y = . sinn (46) 


The folloving relations apply to the large triangle (ABC) in Figure 6. 


Up = |Wp | sin (a + Y )/sin UÇ (47) 
2 2 2 

U,- G^- | Wa | 226 [W, | cos (a +7) (48) 

Iw | = uc بس‎ G cos U - (49) 
h h h o ” 


From equations (31), (45), (46), (47), (48), and (49), the following 
equations describing the surface boundary layer characteristics were 


obtained 
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sinl a + y ) 
T للا‎ = 一 
PO (s/q) - cos ( a + Y) 8 


k Gq sin (a + y) 


U; = s sin u. [ In BRKT U (R/L)] (51) 


By including (23) and (15), Cardone obtained a system of equations that 
could be solved for U, and "بب‎ from the input parameters G, f, 

ӘС | | 
1 ٣ oe); ee and n . These equations were calculated from a 
knowledge of sea level pressure, surface air temperature, sea surface 
temperature, and latitude. 


.. Some Characteristics of the Marine Surface Boundary Layer 





A primary effort of Cardone's work was to determine how 
well the planetary boundary layer model, using only large-scale synoptic 
parameters as input, specified the surface boundary layer wind distribu- 
tion over the sea surface. Cardone considered several studies relating 
the geostrophic wind ratio U,/G and the air-sea temperature difference. 
The results of an extensive study by Carstensen (1967) at Fleet Numerical 
Weather Central (FNWC) appear in Figure 7 along with earlier results 
due to Bleeker. The ratio of surface wind (V) to geostrophic wind (G) 
is tabulated for each degree of the reported air-sea temperature differ- 
ence, and the median value for each interval is plotted as an x in 
Figure 7. The validity of the results was limited to the air-sea tempera- 
ture difference range -4 to +1ºC. In order to include observed winds 
in the surface pressure analysis, FNWC uses a constant value for V/G 
of 0.78 and a turning angle of 15º toward lower pressure. This selection 


could produce significant error. 
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Fig. 7. Reported wind-geostrophic vind ratio as a 
function of air-sea temperature difference 
(Carstensen, 1967) 
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Fig. 8. Two layer model predictions of the dependence 
of the ratio of the 19.5 meter wind speed to 
the geostrophic wind speed upon air-sea 
temperature difference and thermal advection 
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Cardone's model can be used to compute a wind speed at 
19.5 meters. Figure 8 illustrates the variation of the ratio of the 19.5 
meter wind and surface geostrophic wind with the air-sea temperature 
difference. The ratio varied less in unstable conditions than in stable 
conditions. 

The effects of baroclinicity also appear in Figure 8 and 


ee 


these effects increase as stability decreases, with the ratio, Uig 5 
increasing with cold advection and decreasing with warm advection. 
Mendenhall (1967) also studied the influence of baroclinicity and sta- 
bility on the veering of the wind in the planetary layer. He concluded 
that over mid-latitude oceans the diurnal variation of frictional veering 
of the wind with height was negligible, but that lapse rate and baro- 
clinicity were important. Results obtained by Cardone, using the plan- 
etary model, showed qualitative agreement with Mendenhall's conclu- 
sions. Cardone's results also suggested that baroclinicity was more 
significant than stability in determining the veering of the wind in the 
planetary boundary layer. For both cases in situations of strong, cold 
advection the actual wind may back with height, especially in conjunc- 
tion with unstable stratification. 

Additional results on the influence of baroclinicity came 
from a study by Clarke (1970), who found a tendency for the boundary 
layer wind to back instead of veer with height in die conditions. 
He stated that this could have been due to baroclinicity or to accelera- 


OT 
tion. He noted that it could have occurred if = : was systematically 
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Ux 
negative in convective conditions up to a level of about 0.2 ES 


instead of being positive as in the normal case. This latter interpreta- 
tion was based on the observation that the wind appeared to be sub- 
geostrophic in the convective layer and reached geostrophic values only 


after a rapid increase through the overlying inversion. 
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HI. COMPUTATIONS AND NUMERICAL PROCEDURES 


A. EXTERNAL WIND FIELD SPECIFICATIONS 

All computations were made with data defined on a subset of the 
FNWC 63X63 grid which covered the North Pacific Ocean (Figure 9). 
Data fields included sea level pressure, sea surface temperature, and 
surface air temperature at each of 315 grid points spaced 381 km apart 
at 60°N. The necessary map factor was included for a 1:30 million 
polar stereographic projection. 

The coriolis term appears as an external parameter and therefore 
poses a limitation on the lowest geographical latitude where the geo- 
strophic relation would be valid. It has already been established that 
the dynamic boundary layer was formed by the surface stress and the 
coriolis acceleration. Pushistov (1970) concluded that manifestations 
of the nonlinearity in the dynamics of the boundary layer and the appear- 
ance of a perceptible ageostrophic component of velocity in the free 
atmosphere were to be expected near latitudes of 7-10°. Sheppard 
(1970) suggested 20°N as a limiting latitude for application to boundary 
layer models. A limiting value of 15° latitude was used in this program 
as the limit on geostrophic flow. For latitudes below 15°, the coriolis 
parameter would be kept constant. 

Ec Geostrophic Wind Computations 

Ihe geostrophic wind components were computed from the 
surface pressure field (P) as follows 
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using a centered explicit finite difference technique. The density was 


determined from the equation of state 





)54( وج - يم 


where Ra is the gas constant for dry air, and 8 the surface air tempera- 
ture in degrees Kelvin. The magnitude of the geostrophic wind was 


computed at each grid point from 


2 
ici . . (55) 
g g 
2. Correction for Curvature 


A gradient wind was calculated to correct for the effect of 
curvature on the surface geostrophic wind. Two basic principles were 
required for this correction. The first was that following an air parcel 


along an isobar, the total derivative was zero or 
dP v) m 1: (56) 


The second requirement was the mathematical statement of curvature 


икисини 
Ки" ۶۹۶ 


E 
E (57) 


where y' and y" were the first and second derivatives with respect to x, 


and R was the radius of curvature. By expanding the total derivative of 
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P, obtaining y' and y", and rearranging, the following result was ob- 


tained for the radius of curvature 


5 
ey) ax? | /2 


9P/3‏ 7۳ + +ؤ ۱ ۹ھ" 





2 2 2 (58) 
ar |, ap _ BP max | ӘР ӘгР 
Ox Ox oy 3 y 2 OP/dy öy dx? 
Next the gradient wind was computed from 
y2 
Һә кә (59) 
İR | 


Solving for V_ and rearranging led to a workable form of the equation 


G 
G 


G ` 1⁄2 + A 1/4 + G/R£ ` 0 


For this case the magnitude of the term G/Rf determined how much curva- 


V 


ture would be allowed. The maximum values were chosen to be -3/16 
for anticyclonic curvature and +3/4 for cyclonic curvature which yielded 
a range of VG = 1.33G for anticyclonic curvature and Ve” 0.67G for 
cyclonic curvature. This agreed with generally accepted synoptic limits. 

Subroutine GRADWD (Appendix A) makes the necessary com- 
putations and follows a similar scheme utilized by FNWC (Kaitala, 1971 
unpublished). Centered finite difference schemes were used to calculate 
the derivatives. 

3. Thermal Wind Correction 
The nondimensional thermal wind components were deter- 


mined from 


nə oe o (61) 





25 8 


7 DO 
T t^g, ۵ X 1 ١ (62) 

and the total thermal wind was given by 
كن )و - ع‎ + 05 . (63) 


For this model, r had a limiting value of 100. The geostrophic and 
thermal wind angles were calculated from the previously computed wind 
components. This procedure was repeated for each grid point and was 


used as input data to the boundary layer model. 


Be PROCEDURES IN SOLVING THE BOUNDARY LAYER EQUATIONS 
]. Initial Conditions 
a. Geostrophic Wind Speed Less Than Or Equal to 5.0 
meters/second: 

This case was associated with the aerodynamically 
smooth flow region described earlier. For this case no corrections were 
made for stability and a constant drag coefficient of 0.022 was assumed. 
This resulted in a constant inflow angle of 15?, The 19.5 meter wind 
was determined from 


Une: = 0.7G (64) 


The choice of the drag coefficient was not arbitrary, but depended on the 
relationship between C, and سلا‎ (see Figure 2). The value of Cio was 


a minimum at 6 meters/second under neutral conditions. The value of 


C19 5 was extrapolated from this curve. 
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b. Geostrophic Wind Speed Greater Than 5.0 meters/ 
second and Absolute Value of Air-Sea Temperature 

Difference Less Than or Equal to 1°C: 

Initially the drag coefficient was C, = 0.0245, with 
the thermal wind (r) and the angle between the geostrophic and thermal 
wind (n) set equal to zero. For this case L,, L and U (h/L') were 
set equal to zero and 0 (h/L')21.0. A solution was then obtained using 
the modified baroclinic model of Blackadar described below. 

E Geostrophic Wind Speed Greater Than 5.0 meters/ 
second and Absolute Value of Air-Sea Temperature 

Difference Greater than 1? C: 

The drag coefficient С̧ = 0.0245 is used to estimate 
the initial U, and r and N are set equal to zero. A - p was restricted 
to the range -15 to +4.0°C. The values of L,, L', W (h/L'), and 
B (h/L') were computed from z 0 š 7 and U, and a solution was 
obtained using the modified Blackadar model described below. 

L, was computed from (34) and L' from (15) using as 


a first approximation 


(007) = 10: m (48/2) (65) 


U (h/L') was computed from SUBROUTINE PSI (Appendix À). If L, was 
greater than zero (stable case), then 


W(h/L') = 0.7L, (66) 
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If L, was less than or equal to zero (unstable case), a value of 


Richardson number and 7: (h/L') had to be computed. This was accom- 


plished in SUBROUTINE SHR (Appendix A). Then W (h/L') was calculated 


trom (13). 
pu (h/L') was computed using SUBROUTINE SHR 
(Panofsky, 1968; Cardone, 1969). If L, was greater than zero (stable 


case), then 
9 (h/L') = 1+7L, (67) 


If L, was less than or equal to zero (unstable case), then a new 


Richardson number was determined from 


1/4 


Ri, - L, (1- 18R ) (68) 


ew * lolg 


where initially Rio = L This iteration was continued until the 


ld 


inequality 


ки .. (69) 





was satisfied, where € — 0.01 for this model. Then o (h/L') was 


found from 


"+00872 
0, VE) ИЛ тәр 


(70) 
İnev ) 1/4 


2. The Modified Blackadar Baroclinic Boundary Layer Model 


The basis of this model was an iteration procedure which 


included the effects of stability anc baroclinicity on the computed 


wind. The procedure was free of computational instability and converged 


rapidly for all reasonable choices of initial U, . First values of zo, 


il 





72725 “€” €, Y, Uo, and U, were obtained using respec- 
new 


tively equations (23), (35), (38), (43), (44), (45), (46), (50), and (51). 


The iterative procedure was repeated until the inequality 


Ver - Ux < € (7 1) 





was satisfied where « was 0.05. The iteration was repeated, setting 


U,-U, , until the inequality was again satisfied. If L, was not 
new 
zero, then new values of W, Zr and L' ow were computed where the 
argument of W was z,/L'. If 
1 = ! > 
Һу. L 2 5 (72) 








then L' was set equal to Lona The model was repeated from where 


ə 
L, vvas first calculated until the inequality vvas satisfied. After com- 
pleting the iterations for the case where the thermal wind (r) and n 

were equal to zero, the computed values ofr and N were substituted 
and the model was again repeated until all criteria were satisfied. If 


L'#0, then @(z/L') was computed as £(19.5/L') ; otherwise 0 (19.5/L') 


was set equal to zero. Finally the wind at 19.5 meters was determined 


from 
U 19.5 19.5 
U ہے‎ İn له‎ - y(+= 73 
105 E | Zo 7 ) (73) 
Values of 8 -0 , G, z,, z/L', W,, U,, Ujg 5. 109R? UC A 


О ال‎ and R, were then available for each grid point. 


42 





C. COMPUTATION OF INDEX OF REFRACTION CONSTANT 

In this section possible relationships are considered between the 
refractive index structure constant and quantities dependent only on the 
wind speed and air-sea temperature difference. Specifically those 
quantities that can be predicted from the model considered in this study 
are examined. 

The propagation of light waves and radio waves is affected pri- 
marily by air-density inhomogeneities and therefore by temperature fluc- 
tuations. One effect of these fluctuations is scintillation or rapid changes 
in intensity at a point. These are distortions caused by variations in the 
speed of propagation through turbulent elements having different densi- 
ties, and hence different indexes of refraction. The refractive index 
structure constant (c 2) is the primary parameter characterizing the turbu- 
lent fluctuations of the atmospheric refractive index for scales from a 
few millimeters to several meters. If the quantity c is known, cal- 
culations can be made on scintillation, beam spreading, and beam wander, 
ail of which are important in describing light wave propagation. 

A relation necessary in the formulation is the one-dimensional 
spectrum for temperature fluctuations in the inertial subrange described 
by Corrsin (1951) 


=) 7 
B, (K)= Be /3 4/2 (75) 


where € is the rate of dissipation of turbulent kinetic energy, Nisa 


measure of the "smearing" of temperature inhomogeneities, Bis an 
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empirical constant (20.81, Boston, 1970), and K is the streamwise 
component of wave number. The structure constant for temperature fluc- 


tuations (Cz) IS defined ds 
-1 
2 
C. = B € /3 N (76) 


and is related to the index of refraction as follows 


-6 2 
2 -79 x 10 P 2 
С̧ = (一 一 一 一 ےم‎ 
E ( x ) T (77) 
a 


Another expression which has to be considered is the nondimen- 
sional balance equation for turbulent kinetic energy, which, for station- 


ary, horizontally homogeneous flow is (Busch and Panofsky, 1968) 
TE = = = 
p, E ир dp = 0. )78( 


9 is the nondimensional wind shear (8), Bo represents the diver- 


gence terms, and 
£ kz 
50 nz (79) 
Ur 


The following relations arise from the definitions in the non- 


dimensional turbulent kinetic energy balance and definition of N 


3 
یلا 


KZ 





£ = f (80) 


£ 

_2 
HO ٢ 28 
IZ | Ен اچ‎ 


Substitution of these two expressions, plus the definition of Km (29) 


)81( ل 


into the definition of p yields (Panofsky, 1968) 


-1/3 KH ] 4/3 ad 


2 zın 98 
Cn = B f. CK. ) 8, (kz) (5 11 (82) 
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Equation (82) is not a suitable expression because of the lack of 
information on ¢ . However if the effect of the divergence terms 


is neglected, pe can be expressed in terms of B and z/L as 


E. = Д - Z/L. | (83) 


u 

It is noted that with wind-wave interaction, the assumption that the 

divergence terms can be neglected is probably not valid, but the possible 

effect this assumption has on predictions of e has not been studied. 
Substitution of (83) into (82) yields 


x. 4/3 


2 2 Ser 
С̧у = B Le 2/1] K 3 ع()‎ 2( ки, (84) 


m 


This final equation for Ce lends itself to the boundary layer model be- 
cause all the parameters can be computed within the model. Thus, the 
model can be used to obtain a horizontal distribution of É over water 
for various wind and stability conditions. The average stability 90/ 7 
is determined from (14). 

The ratio of the exchange coefficients KK. is set equal to 1.0 
in Cardone's formulation. However it has been shown that Кр апа K. 
are not usually identical except for small values of 90/ dz for which 
the turbulent motion is not significantly influenced by buoyancy forces. 
Businger et al. (1971) suggested the following form for the ratio of eddy 


diffusions 
1/, 
262 
a rıb Каи 


pL) 85 
h (1- 15 و‎ 0 
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Here a, = 1.35 for neutral conditions, which agrees closely with labora- 
tory measurements (cf., Hinze, 1959). A value of 1.35 rather than 1.0 
for a, under neutral conditions arises because the von Kármán constant 


was equal to 0.35 instead of the usual 0.40 (Businger et al., 1971). 


D, WAVE HEIGHT ANALYSIS 

Wave analysis is another application for the boundary layer model. 
For this purpose, a simplified version of the FNWC sea-swell model 
(Hubert and Mendenhall, 1970) was utilized. The u and v components 
of the 19.5 meter wind determined from the boundary laver model were 
used as both current and history input in this analvsis. 

The FNWC sea-swell model computes a weighted wind from the 


wind components as follows: 


u 0.3 F 0.7u, 


i 


“ии, O (86) 


t 
where t-12 is the 12-hour old history wind speed and t is the current 
wind speed. These component values are saved as history for the next 
analvsis run, and then the total wind (V) is computed. The wave heights 


(WA) in feet are then computed from 


2 
V 
NES 7 Бә 


Although this is a very simple empirical method for determining wave 


5 
V 
= Ga | + 1.0 (87) 


heights, the basic method has been used by FNWC since the mid-1960's 


and this particular model has been used operationally since August 1970. 
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In addition, the FNWC model computed corrections for fetch and 
ice, and performed some smoothing of the resulting height field. It was 


assumed these corrections would be minimal for this application. 
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IV. SPECIFICATION OF EXTERNAL PARAMETERS 


A primary objective of this study was to apply the two-level baro- 
clinic boundary layer model to actual atmospheric conditions and compare 
the computed results with the observed results. The time period 00 GMT 
28 November 69 through 00 GMT 02 December 69, with an analysis inter- 
val of every 12 hours, was chosen for this comparison for several reasons. 
First an intense low pressure system, which influenced most of the central 
North Pacific Ocean, developed during this period with the system reach- 
ing maximum intensity about 00 GMT 30 November 69. The system pro- 
vided input data which allowed the model to predict boundary layer 
parameters under varying atmospheric conditions. Secondly considerable 
synoptic data had been gathered by the author during an earlier study on 
the extreme sea conditions (unpublished) associated with this storm. The 
latter gave a good basis for comparing computed values to observed data. 
Finally, this storm produced waves and swell which caused considerable 
damage to several islands of the Pacific as well as to the spar buoy FLIP. 
Therefore a wave generation model was applied using the predicted sur- 
face winds, and these predicted waves were compared to the observed 
wave conditions. 

The 00 GMT 30 November 69 period was of the most interest due 
to its extreme conditions, extensive data coverage, and it was the prob- 
able period when the destructive waves were being generated. The various 


corrections described above were incorporated into the model for this 
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time period, and synoptic data were plotted and analyzed for surface 
pressure, air-sea temperature difference, surface wind speed and wave 
height. 

Surface layer parameters were computed including various combina- 
tions of stability, curvature, and thermal wind to yield information on 


the relative importance of each parameter in the formulation. The follow- 


ing notation will be used in describing the results: 


NN - non-neutral, correction for stability 
N = neutral = no stability correction 

NC - no curvature correction applied 

C = curvature correction applied 

NT - no thermal wind correction applied 
T - thermal wind correction applied 


It was recognized at the onset of the study that the specification 
of the external fields ( A9 and G) had to be reasonably accurate in 
order to evaluate the model on the basis of observed wind fields. For 
example, if the observed surface winds were quite different than those 
predicted by the model, some estimate had to be available on the possi- 
bility that the discrepancies were due to poor specification of the geo- 
stophic wind or the air-sea temperature difference. Therefore in this 
section consideration is given to the relevant aspects of representing 
those fields which make up the external parameters. 

The charts in the following figures correspond to a 1:30 million 


polar stereographic projection with the grid as depicted in Figure 9. 
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A, FNWC ANALYSIS 

Sea level pressure and air-sea temperature difference fields for 
00 GMT 30 November appear in Figure 10. A 967 mb low pressure 
center was located near 45°N 180° with considerable thermal instability 
both ahead of and behind the system. A significant feature of this sys- 
tem was the intense isobaric gradient in the western quadrant. 

The first consideration was the comparison of the computed surface 
geostrophic wind with the m cs gradient wind for this period. The 
geostrophic and gradient wind fields computed from the surface pressure 
fields in Figure 10 appear in Figures 11 and 12 (units of meters/second). 
Figure 13 depicts the influence of the curvature correction on the geo- 
strophic wind computed as a percent difference from geostrophic wind 
minus gradient wind then divided by geostrophic wind. 

The two maximum wind centers west and north of the Low are based 
on data at one grid point and result from the influence of anticyclonic 
curvature, strong isobaric gradients, and finite difference computations. 
These features do not appear to be realistic. However, the wind speed 
distribution depicted by them is good except that the magnitude of the 
wind speed may be a little too large. It was expected that the maximum 
change would occur in an area of maximum curvature and maximum pres- 
sure gradient from the relation in (60). Generally speaking, the change 
in wind speed due to the curvature correction is 10 to 20% with a maxi- 


mum change of 34% imposed by the limits placed on the term G/Rf in (60). 


50 








٥ا٥1‎ 018106 285-182 pue (pTlOS) eunsseud eoeJuang 


uo peseq 69 AON 


a * \ 
59, 





BEP IANA 
ОС IND OO 4ОЈ 9ртәтЈ (рәддор) 


` 


OLL 





"Ol 


`x ° 
N 


PO 
Sc, 
۹ 


STA- 











LU. اا‎ 
= 2" 7. 








( 


8STsÁ[vue poo 


Ol əng) eunssoeud 
ds puTm 51110404505 


x 


“STA. 





52 








(Ot eunstwy) eunsseud 
soejans uo peseq 8٢٥9۸٢608 5 ритм 10٥179٩٥ “2! 


"TJ 


.. 





sik 





-и). 
ok” 
7 


рд 5. 
: | Y ки” ( 
o Oss 1 
175  هہب‎ ٦ 


. 
t 


G, 


.9 . 


4ueTpuddS-oTudou48s 





( 


O 


Cl eundTH-=l| eundTH) ертәтј ритм 
әј TOT PTOTI SdUSISIJTP ۹1393 dəq "Çi 












“204. 












| " = m | u 1 
p na 0 = PE 


| ” иә = yə — ə 
73 ۰ | : o a 2 E x PS 1 











B. SUBJECTIVE ANALYSIS 

The results from the hand analyzed surface chart were also used 
in comparisons of results. Figure 14 illustrates the surface pressure 
and air-sea temperature difference fields from the subjective analysis 
which are analogous to Figure 10. The major variations between Figure 
10 and 14 appear to be due to the difference in the number of synoptic 
reports used as wellas the analysis scheme. Theoretically, Figure 14 
should be more representative since more data were included, but it was 
a subjective analysis. Figures 15 and 16 depict the geostrophic and 
gradient surface winds computed from the subjectively analyzed pressure 


field. 
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V. RESULTS FROM EXPERIMENTS WITH THE MODEL 


A. COMPARISONS IN TERMS OF SURFACE ROSSBY NUMBER AND 

STABILITY CATEGORIES 

Near surface parameters were computed for the 19.5 meter leve1. 
"This level represents an average height of ship anemometers (Bunting, 
1968). Theoretically the computed 19.5 meter level winds could be 
compared with the observed wind fields and this would vield some 
information on the accuracy of the model in predicting surface winds. 

In order to verify the computations and also to examine the range 
of conditions provided by the data, a comparison was made with Cardone's 
results which appear in Figures 4 and 5. Results for this comparison 
were obtained from computations in which stability was included but not 
baroclinicity. Results from all nine time periods appear in Figures 17 
and 18 where values of the logarithm of the surface Rossby number (24) 
range from 7.76 to 9.45. The U,/G versus log RO results (Figure 17) 
agree with Cardone's results (Figure 4) except with regards to the inter- 
valbetween L, isolines. In Figure 17 this interval (spacing) is larger. 
With respect to > versus log Ro (Figure 18), the unstable results 
(negative Ly) agree with those in Figure 5 except, again, with regards 
to the interval between L, isolines. The stable results (positive L,), 
however, deviate considerably from Cardone's results with the veering 


angle decreasing as L, increases. The results in Figure 18 would 
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indicate that the veering angle decreases with both increasing stability 
and increasing instability with a discontinuity occurring near neutral 
conditions. 

Ihe unusual, and unexpected, features in the results just de- 
scribed were examined in considerable detail. The same results were 
noted by Pafias (unpublished study, 1971) in an investigation with 
Cardone's program. These features in the results are perhaps related 
to the drag coefficient representation selected by Cardone and also be- 
cause different iteration limits were used in this study than those used 
by Cardone. 

The range of data with respect to U,/G and 9 - 9 was also 
examined and the results appear in Figure 19. The computations for 
these results were made on the basis of a barotropic atmosphere with no 
correction for isobaric curvature. These computations differ from those 
used by Cardone for Figure 8 in that a variable coriolis parameter was 
used. Isolines have been drawn for geostrophic wind speeds ranging 
from 15 to 75 knots. For those wind speeds considered by Cardone in 
Figure 19, and also for wind speeds above 20 knots considered by 
Pafias, these results are in agreement. The change in shape of the 
curves for winds under 20 knots was not described by Cardone but is 
probably due to a combination of the functional forms for the stability 


dependence and the drag coefficient as chosen by him. 


62 





0.9 


0.8 


0.7 


Och 


15 


-7.0 -5 ә0 - 3,0 
Ва ~ Os 


-] әО 


50 


70 


1.0 


Fig. 19. Geostrophic drag coefficient versus air-sea 


temperature difference for range of geostrophic wind 
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B. EXPERIMENTS WITH CHANGING THE PHYSICS IN THE TWO-LAYER 
MODEL AND CHANGING THE SPECIFICATION OF THE SURFACE 
GEOSTROPHIC WIND 
In this section Ug 5 results which were obtained from different 

computations with the two-layer model are compared. These computations 

were varied by including or not including specific parameters or correc- 
tions. The variations considered and compared were described earlier 
and are 1) N/NC/NT, 2) NN/NC/NT, 3) NN/NC/T, 4) NN/C/T, and 

SENN/NC/T with k=0.35. 

These variations will be referred to as Cases 1 through 5. Itis 
noted that if the geostrophic wind speed was less than or equal to 5.0 
meters/second, stability and baroclinicity were not taken into account 
in the iteration schemes. In regions with this wind speed range, espe- 
cially in the subtropics, differences will not arise except where isobaric 
curvature exists. 

Per cent difference (1-2) will be used to compare the two cases 


where (1-2) corresponds to the following computation at each grid point: 


case l result - case 2 result 
case 1 result 


X 100 

Although large differences between cases were unexpected, small cor- 

rections at low wind speeds could produce large per cent changes. 
The surface wind fields corresponding to each of the five cases 


are presented in the following paragraphs along with selected comparisons 


between pairs of cases. 
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Case 1 (N/NC/NT) 


This is the simplest case with respect to the possible physics 
which are included in the model. The external wind field is defined for 
a neutral, barotropic atmosphere with no curvature correction. The wind 


field at 19.5 meters for this case appears in Figure 20. 


Case 2 (NN/NC/NT) 


The effects of stability are now included and the 19.5 meter wind 
Held for this case appears in Figure 21. The per cent difference (2-1) 
field is shown in Figure 22. Positive values correspond to increases in 


the computed U field due to including the stability effect. This in- 


1923 
crease is quite evident on the western side of the cyclone. The patterns 
are not expected to be coincident with the air-water temperature isolines 
in Figure 10 because the stability effect is dependent on both the temper- 
ature gradient and the wind speed. The 20% contours near the eastern 
border are attributed to small corrections at low wind speeds as dis- 
cussed previously. 

Case 3 (NN/NC/T) 

The effects of including the baroclinicity or thermal wind are now 
included and the 19.5 meter wind field for this case appears in Figure 
23. The per cent difference (2-3) field is seen in Figure 24. Negative 
values correspond to an increase in the computed 19.5 meter wind due 
to including the thermal wind in the model. In this difference field, only 
20 grid points (6%) appeared to be significantly affected by including the 
thermal wind. These grid points generally lie in an area of warm air 


advection. as 
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Case 4 (NN/C/T) 

In addition to all specifications possible within the physics of the 
model, the correction for curvature is applied to the external wind field. 
The 19.5 meter wind analysis for this field is shown in Figure 25 and 
the per cent difference field (3-4) appears in Figure 26. Negative values 
correspond to larger 19.5 meter winds due to including the curvature cor- 
rection. Strong wind maximums, noted in comparing the external wind 
fields with and without this correction, also appear in Figure 26 as the 
predominant feature. In general, the effective change in the 19.5 meter 
wind field specification is larger and covers more area than any of the 


other cases. 


Case 5 (NN/NC/T with k= 0.35) 


This case differs from case 3 only with respect to the change of 
the von Kármán constant from 0.40 to 0.35. The 19.5 meter wind field 
for this case appears in Figure 27 and the per cent difference field (3-5) 
is depicted in Figure 28. Negative values correspond to increased 19.5 
meter winds due to the change in the value of k. The suggested value 
of k=0.35 arose from recent surface layer observational investigations 
by Businger et al. (1971). They noted from their results that using 
k=0.40 in several empirical expressions led to estimates of U, values 
which were 15% too high when compared with direct measurements of 
the momentum transfer. The fields in Figure 28 indicate that a general 
increase of 12 to 17% occurs in the computed 19.5 meter wind field if 


the lower value of k is used. This increase agrees with the U, 
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discrepancies noted by Businger et al.. Because the von Kármán 
constant appears in several expressions in the boundary layer model, 
the importance of its change in the two-layer model was not known. 
These results indicate that the difference between k2 0.35 and k2 0.40 
is important, 

In summary, the five cases described above indicate which terms 
were most influential in the model. It should be remembered that most 
of the subtropical points were not influenced by the model since the 
wind speeds were generally less than 5.0 meters/second. Of all the 
“° modifications to the mode], the curvature correction to the external 
field produced the largest variation over the entire grid. Changing k 
from 0.40 to 0.35 accounted for the second most significant change. 
Differences in the surface wind fields were evident, however, when 
specifications in the model were made. Of the two possible specifica- 
tions, stability appeared to be more significant than the thermal wind 
with respect to area influenced and per cent difference changes. 

Although the influence of the thermal wind in this model appears 
to be small with respect to the relative magnitude of ine computed wind, 
the influence it has on the direction of the wind appears to be signifi- 
cant. The difference field, in degrees, of the wind direction for 
NN/NC/T - NN/NC/NT appears in Figure 29. This field represents the 
difference in wind direction between a case with and a case without the 
thermal wind specification. The region of significant differences does 


not coincide with that region where the thermal wind effect was large 
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withrespectto wind speed. Positive values in Figure 29 correspond to 
a veering wind and negative values correspond to a backing wind. It is 
noted that the direction change due to baroclinicity is primarily deter- 


mined by the angle between the thermal and geostrophic winds. 


C. COMPARISONS OF COMPUTED WINDS WITH OBSERVED WINDS 

These comparisons will be described for the case NN/C/T with 
k=0.35. The 19.5 meter wind fields computed from the other cases were 
compared to the subjectively analyzed observed surface wind field 
(Figure 30) and in all cases the computed winds underestimate the ob- 
served winds. This was also the case when the above corrections were 
included, but yielded the best results with respect to the magnitude and 
patterns of the wind speeds. 

Because the previous comparisons revealed the importance of the 
specifications of the external wind field with respect to curvature, it 
was considered necessary in this comparison to examine the importance 
of the initial analysis (FNWC or subjective). The 19.5 meter wind field 
computed with the above specifications and with the FNWC analysis is 
shown in Figure 31 and the wind field based on subjective analysis is 
seen in Figure 32. The observed surface wind field which appears in 
Figure 30 was analyzed using only those observations reported above 
10 meters/second. In Figures 31 and 32, the maximum contour, 30 
meters/second, agrees with the observed maximum wind speed, but the 


objective analysis does not indicate the 25 meter/second maximum east 


of the low. 
78 








pezKieus Á[eAT[39e[qnuns pue م2420‎ PO990Td 
srsäTueus (‚955 sm) pəəds pu[A eoveJIJns 


E 5 ۳:٠٢ NÓ 77 


T «c 


. 


UO 655 
рәлләвао *OC *3T4 . 


777: 


ə. )5ہ‎ 


70 


* 
Y 一‏ — €( ہو جو بے ہے په وت اه 


一 一 一 


1 


..........3....—-.......üu. 


7 





سس ج سر ә‏ 


° 
y 
. 


possa “pəpnTouT suio4 110174207370205 9 0 : 

RUE UYTM 8 Siu жҝ Pen Дәфә 2٦ سس‎ ATE ° 
85 "ء۰۰‎ ıı “7 207077. 5٤ 
2035 | J Neer 


Bep OMNA UO 








80 





098 1 


0 


6 





0 


OX 


9180 5 0ہ‎ uo 
02980017 еидәд uoTI1BOTITOOCAS 5118 106 2 


*u) sys4Teue pasds puja Jeqau G° 61 "ət "STA‏ 2998¡( ۸پب۷پ ھ۶ 


2m 5 se 





ә” | 


81 





A comparison of computed 19.5 meter winds from the two analyses 
(Figures 31 and 32) and the external wind fields (Figures 11 and 15) 
shows that differences in the external wind fields are 0 ии 
carried through the model. The model, however, was able to produce 
19.5 meter wind fields in both cases which were in good agreement with 


the observed wind field, at least for wind speeds over 10 meters/second. 


D, RESULTS FROM WAVE HEIGHT COMPUTATIONS 

Wave height fields were computed using the 19.5 meter wind 
fields obtained from the boundary layer model and the simplified wave 
height program developed from the FNWC operational sea-swell model. 
The program was run starting at 00 GMT 28 Nov in order to build history 
for this analysis. The 19.5 meter wind computations were obtained 
using all specifications and corrections described previously, and they 
were also made on both the FNWC and subjectively analyzed fields. The 
computed wave field from the FNWC data is seen in Figure 33 and the 
wave field from the subjective analysis appears in Figure 34. The ob- 
served wave field obtained from synoptic observations is depicted in 
Figure 35. It is noted that wave height is probably the most difficult 
parameter to observe and therefore a representative analysis of the ob- 
served wave field is difficult to achieve. 

The wave field obtained from the subjective analysis, Figure 34, 
appears to be in good agreement with the observed heights with maxi- 


mum computed heights slightly less than the maximum observed heights. 
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Results computed from the FNWC analysis, Figure 33, appear to have 
similar agreement with respect to the maximum height but poorer agree- 
ment with respect to the overall pattern of the wave field. 

In general these results indicate that the two-layer model can 
produce surface wind fields which can be used for wave height analysis, 
even from a very simple wave model. It is noted that in the subtropics 
the wave heights computed from both sets of analysis are considerably 
below the observed wave heights. This is probably due to low surface 
geostrophic wind speeds in this area and the fact that no corrections are 
made for stability or baroclinicity at low wind speeds. This could pos- 
sibly be corrected in the model by selecting a higher initial value for 


the drag coefficient. 


E RESULTS FROM INDEX OF REFRACTION COMPUTATIONS 

A final application of the boundary layer model was the computa- 
tion of the index of refraction structure constant (Cİ) using (84) mne 
ratio of the exchange coefficlents was determined from (85) and the 
input data for this equation were obtained by including the stability 
effect, the therma] wind effect, the correction for curvature, and using 
k=0.35. All parameters needed to compute on were determined in the 
model. 

The results from these computations appear in Figure 36. Values 
әр с | range from 0.0 in regions of near neutral stratification to a maxi- 


1 


mum of 0.181 in the region of maximum wind speeds. Comparisons of 
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these values with those measured by Portman et al. (1968) indicate that 
they are comparable in regions of low wind speed with = - ә. less than 


] degree. Although no definite relation between ” and wind speed or 


E - ә. appears in Figure 36, larger values of C x appear to be associ- 


T 
ated with the greater wind speeds and larger instability (9. - 9 ranging 
from -1 to -6?C). 

The results indicate that a horizontal distribution of E for a 
given height, in this case 19.5 meters, can be computed with the model 
with some degree of confidence in the range of values obtained. Com- 


parisons should be made between results computed by this model and 


observations of C. to obtain better verification. 
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VI. CONCLUSIONS 


The results of this study have indicated that a marine boundary 
layer model can be applied to synoptic scale data to obtain various sur- 
face parameters. These parameters include friction velocity 07 
surface roughness (25) , surface Rossby number (Ro? , 19.5 meter wind 


(U ), veering angle (Wo), and stability length (z/L). 


0+025 

First, the results indicate that the curvature correction produced 
the greatest change to the wind analysis. However, stability and modi- 
fication of the von Kármán constant produced significant changes in the 
field, while changes due to the thermal wind were confined primarily to 
changes in wind direction. 

Second, the 19.5 meter wind was computed including various 
physical specifications and corrections and compared with synoptic 
observations. The results indicate that the computed wind was lower 
than the observed wind in all areas and for all cases, but the general 
patterns agreed well with observations. 

Finally, it appears that Cardone's planetary boundary model pro- 
vides a reasonable framework to specify those surface layer parameters 
required for spectral wave forecasting and index of refraction calculations. 

Further studies are needed to more completely verify and improve 
the model. The specification of the drag coefficient is probably the 


aspect in this model which needs the most attention. The drag coefficient 
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has received considerable examination in observational investigations 

and the results indicate some disagreement among investigators. A 
suggestion of this study is to increase the initial value of C, to increase 
the computed surface wind speeds especially in areas of low wind speed. 
Application of this model to a region of reliable synoptic data such as 

. during the BOMEX experiment in 1969 would reduce some of the errors 
due to data input and also test the physics of the model in the tropical 
atmosphere. Finally, observations of the index of refraction over a wide 
range of wind speed and stability would allow reliable comparisons to be 


made of the structure constant (Cr). 
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